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We present a systematic theoretical study of the BCS-BEC crossover in two-dimensional Fermi gases with 
Rashba spin-orbit coupling (SOC). By solving the exact two-body problem in the presence of an attractive short- 
range interaction we show that the SOC enhances the formation of the bound state: the binding energy E B and 
effective mass m 3 of the bound state grows along with the increase of the SOC. For the many-body problem, 
even at weak attraction, a dilute Fermi gas can evolve from a BCS superfluid state to a Bose condensation of 
molecules when the SOC becomes comparable to the Fermi momentum. The ground-state properties and the 
Berezinskii-Kosterlitz-Thouless (BKT) transition temperature are studied, and analytical results are obtained in 
various limits. For large SOC, the BKT transition temperature recovers that for a Bose gas with an effective mass 
niB- We find that the condensate and superfluid densities have distinct behaviors in the presence of SOC: the 
condensate density is generally enhanced by the SOC due to the increase of the molecule binding, the superfluid 
density is suppressed because of the non-trivial molecule effective mass m^ . 
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It has been widely believed for a long time that a smooth 
crossover from Bardeen-Cooper-Schrieffer (BCS) superfluid- 
ity to Bose-Einstein condensation (BEQ of molecules could 
be realized in an attractive Fermi gas fljj— 13|] . This BCS-BEC 
crossover phenomenon has been successfully demonstrated 
in ultracold fermionic atoms by means of the Feshbach res- 
onance |4]. Some recent experimental efforts in generating 
synthetic non-Abelian gauge field has opened up the oppor- 
tunity to study the spin-orbit coupling (SOC) effect in cold 
atomic gases |5J]. For fermionic atoms jgj, it provides an al- 
ternative way to study the BCS-BEC crossover |01 according 
to the theoretical observation that novel bound states in three 
dimensions can be induced by a non-Abelian gauge field even 
though the attraction is weak |8l |9[] . 

Recently, the anisotropic superfluidity in 3D Fermi gases 
with Rashba SOC has been intensively studied | nJ-12]. Two- 
dimensional (2D) fermionic systems with Rashba SOC is 
more interesting for condensed matter systems II 1311 and topo- 
logical quantum computation fl4ll . By applying a large Zee- 
man splitting, a non-Abelian topologically superconducting 
phase and Majorana fermionic modes can emerge in spin- 
orbit coupled 2D systems 11411 . In the absence of SOC, 
the BCS-BEC crossover and Berezinskii-Kosterlitz-Thouless 
(BKT) transition temperature in 2D attractive fermionic sys- 
tems were investigated long ago 1 15 . 1(1 (see lfl7ll for a re- 
view), which provide a possible mechanism for pseudogap 
formation in high-temperature superconductors 11811 . 

In this Letter we present a systematic study of 2D attrac- 
tive Fermi gases in the presence of Rashba SOC. The main 
results are summarized as follows: (i) The SOC enhances the 
difermion bound states in 2D. At large SOC, even for weak 
intrinsic attraction, the many-body ground state is a Bose- 
Einstein condensate of bound molecules. In the presence of 
a harmonic trap, the atom cloud shrinks with increased SOC. 
(ii) The BKT transition temperature is enhanced by the SOC 
at weak attraction, and for large SOC it tends to the critical 
temperature for a gas of molecules with a nontrivial effective 



mass. The SOC effect therefore provides a new mechanism 
for pseudogap formation in 2D fermionic systems, (iii) In 
the presence of SOC, the superfluid ground state exhibits both 
spin-singlet and -triplet pairings, and the triplet one has a non- 
trivial contribution to the condensate density. In general, the 
condensate density is enhanced by the SOC due to the increase 
of the molecule binding. However, the superfluid density has 
entirely different behavior: it is suppressed by the SOC due to 
the increasing molecule effective mass. 

Model and effective potential — A quasi-2D Fermi gas 
can be realized by arranging a one-dimensional optical lattice 
along the axial direction and a weak harmonic trapping poten- 
tial in the radial plane, such that fermions are strongly con- 
fined along the axial direction and form a series of pancake- 
shaped quasi-2D clouds | T^- 21 ]. The strong anisotropy of 
the trapping potentials, namely cj z » a>± where a> : (oj±) is 
the axial (radial) frequency, allows us to use an effective 2D 
Hamiltonian to deal with the radial degrees of freedom. 

The Hamiltonian of a spin- 1/2 attractive Fermi gas with 
Rashba SOC is given by H = J d 2 riff(r) (<H + <H S0 ) if/(r) - 
U j ^r^ T (r)^;(r)^;(r)^ T (r), where if/ = [^ T , ^] T represents 

the two-component fermion fields, 'Ho = — - - hcr z is 
the free single-particle Hamiltonian with // being the chemical 
potential and /; the Zeeman spli tting , and 'K so = —i%A(o~ x dy — 
o~ y d x ) is the Rashba SOC term 12211 . Here cr A v z are the Pauli 
matrices which act on the two-component fermion fields. The 
short range attractive interaction is modeled by a contact cou- 
pling U 112311 . In the following we use the natural units 
% = &b = m = 1. 

In the functional path integral formalism, the partition func- 
tion of the system is Z. — f !Dif/!Diff exp {—S[iff, iff]}, where 

S[if/,tft] = j£ dr[J d 2 rif/d T iff + H (iff, if/)] with the inverse 
temperature f3 = l/T. Introducing the auxiliary complex 
pairing field O(jc) = -Uif/i(x)if/^(x) [x = (t, r)] and applying 
the Hubbard-Stratonovich transformation, we arrive at Z = 
j DyD*¥D<I>D<&* exp {± fdxf dx r ¥(x)G- l (x, x'y¥(x') - 
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U~ l f dx\<&(x)\ 2 }, where <F = |>,iA] T is the Nambu-Gor'kov 
spinor. The inverse single-particle Green function G~'(x, x') 
is given by 



~d T -7/ -74 
-/cr v O*(jt) 



i'cr v O(jc) 
-d T +<H -<H' 



S(x-x'). (1) 



Integrating out the fermion fields, we obtain X = 
J DODO* exp{ - .SefftO, O*]}, where the effective action 
reads S eff [0,0*] = ET 1 J ^|0(x)| 2 - iTrlntG" 1 ^^)]. 

Two-body problem — The exact two-body problem at van- 
ishing density can be studied by considering the Green func- 
tion T{Q) of the fermion pairs, where Q = (/v„,q) with 
v„ = 2rmT (n integer) being the bosonic Matsubara fre- 
quency. In the present formalism, T~ l (Q) can be obtained 
from its coordinate representation defined as r~ (x, xf) = 
(J3V)- l 6 2 S eS [<b,®*]/[6®*(x)6®(x')]\® =l) . For = 0, the 
single-particle Green function reduces to its non-interacting 
form g (K) = diag[g + (K), g_(K)] with g ± (K) = [iu„ + 
(£k _ ho-,) - A(cr x ky + cr v fc v )] _I , where K = (ico n ,k) with 
u>„ = (2n + 1)jtT being the fermionic Matsubara frequency. 
Here ft — - p and e k = k 2 /2. The single-particle spectrum 
generally has two branches: co^ — ft + V/l 2 k 2 + h 2 . 

After the analytical continuation iv„ — > u>+iO + , the real part 
of T~ l (Q) takes the form 



rVq) 



1 - y 

U t-L, 



i -/(«■£)-/(«■£) 



£y,y-±;k 



-co) 



(l +ayT kq ), (2) 



where f(E) = l/(e^ £ + 1) is the Fermi-Dirac distribution func- 
tion, and 7" kq = (A 2 k ■ p + h 2 )/ ^(A 2 k 2 + h 2 )(A 2 p 2 + h 2 ) with 
p = k + q. T l takes the form similar to that of the rela- 
tivistic systems J3, due to the fact that < H so behaves like a 
Dirac Hamiltonian. Since in 2D the bound state forms for 



arbitrarily small attraction M25I1 . the contact coupling U can 
be regularized by the two- body problem at vanishing SOC, 
t/ 1 = 2k(2fk + ^b) 1 OJE LLZD , where e B is the binding en- 
ergy at vanishing SOC. This equation recovers the exponential 
behavior — 2Aexp(-47r/t/) in 2D l26ll . where A » ee is 
an energy cutoff. All physical equations are finally UV con- 
vergent in terms of eb and we set A — > oo in the dilute limit. 




FIG. 1: The binding energy Eg (left, divided by 6b) and the effective 
mass niQ (right, divided by 2m) as functions of t]/£b- 



From now on we consider the case h — 0. The binding 
energy Eb at nonzero SOC is determined by the solution of 
co + 2p = -Eb for F~'(w, q = 0) = 0. From the imaginary 
part of T~ l (Q), the bound state corresponds to the solution 
in the regime -oo < co + 2p < -A 2 and hence Eb > A 2 . 
Completing the momentum integrals analytically, we obtain a 
simple algebraic equation for Eb 12711 . 
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V£b - A 2 



(3) 



The solution can be generally expressed as Eb = £b + 
477/(77/ €b) where 77 = A 2 /2. For 77 «: €b, we have 7^1 and Eb 
is well given by Eb - £b + 2/1 2 . For 77/eB — > 00 , the solution 
approaches very slowly to the asymptotic result Eb - A 2 . In 
general, Eb increases with increased SOC, as shown in FigQ] 
It is straightforward to show that the bound state contains both 
spin singlet and triplet components @]. 

For small nonzero q, the solution for co can be written as 
co + 2p = -Eb + q 2 /(2rriB), where niB is the molecule ef- 
fective mass. Substituting this dispersion into the equation 
r -1 (w,q) = we obtain (27 1 



K-2 



2 m 1 2 Vk - 1 - (k - 2)(| - arctan 

»Jb 2k 2 V/c - 1 + (5 - arctan -^=) 



) 



(4) 



where k = Eb/A 2 . For A — > 0, we obtain the usual result 
ins — > 2m. For A — > 00, we have Eb — » /I 2 and 77?b approaches 
the asymptotic result 4?7t. In general, ms is larger than 2m, as 
shown in FigQ] Together with the result for Eb, we conclude 
that a novel bound state (referred to as rashbon lioll ) forms. 
It would have significant impact on the many-body problem 
discussed in the following. 

Ground state — For the many-body problem, we con- 
sider a homogeneous Fermi gas with fixed fermion density 
n = N/V. For convenience, we define the Fermi momentum 
via n = k 2 /(2n) and Fermi energy by ep = k 2 /2. The ground 
state (T = 0) can be studied in the self-consistent mean-field 
theory, where we replace the pairing field O by its expectation 
value (O) = A. Without loss of generality, we set A to be real. 

The mean-field ground-state energy Q = S e g[A, A] / (J3V) 
can be evaluated as Q = A 2 /U + (l/2)£k(2ft 
E~), where E* = [(ft) 2 + A 2 ] 



2-11/2 



c k 



cle excitation energies with 
ing to the equation that Eb 



the quasiparti- 
ft* = ft ± <l|k|. Accord- 
satisfies, £1 can be evalu- 
ated as Q = Q2d(A,//, €b) + Qa, where Q2d(A,//, 6b) = 
(A 2 /47r){ln[( jFTtf-n)/eB] - 1/2 - fx/( ^Jhs 
formally the ground-state energy for vanishing SOC Il5lll7ll . 
and Q, = -(A/2n) jf dk[ V(& - rj) 2 + A 2 - (ft - 77)] is the 
contribution due to the SOC effect. 

From the explicit form of the ground-state energy, the gap 
and number equations can be expressed as 



[p 2 + A 2 ) 1/2 - p = e B exp [2h Ou/77, A/77)], 

[p 2 + A 2 ) 112 + p = 2e F - 277 [1 - h Qi/n, A/17)] . (5) 
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FIG. 2: (Color-online) The pairing gap A (left, divided by £p) and the 
chemical potential /i (right, divided by e F ) as functions of A/kp. The 
dashed lines represents the analytical results A = ^2E B £p((k) and 
H = —Eb/2 with £b calculated from Eq. (O. 
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respectively. Here the functions I\ and I2 are defined as 
h(a,b) = J Q l dx[(x 2 -l-a) 2 +b 2 ]- 1/2 wdl 2 (a,b) = dx(x 2 - 
1 - a)[(x 2 - 1 - a) 2 + b 2 ]~ 1/2 . I\, h and £l A can be analytically 
evaluated using the elliptic functions. For vanishing SOC, we 
recover the well-known analytical results, A = yfle^ep and 
fi - £f - e B /2 115D . 

Now let us start from weak attraction, e% <sz ep. For suffi- 
ciently small SOC, we have /[ — > and h — * - 1 , and the solu- 
tion is well approximated by A ^ -\f2egep and ep - £b /2 - 
2t], which indicates a BCS superfluid state. For large SOC, we 
expect that p. becomes negative and \fx\ » A. Substituting this 
into the gap equation, we find p -Eb/2, which indicates 
a Bose-Einstein condensate of molecules with binding energy 
Er. Then expanding the number equation in powers of A/\p\ 
and keeping the leading order, we obtain A ^2Eb£f((k), 



where £(*■) = 2k~ 1 (k - l) 3/2 (2 - 1 + § - arctan ; 



K-2 



r 



2 2V«-1' 

This is a transparent formula to show that the pairing gap A in- 
creases with increased SOC, consistent with the perturbative 
approach 112 811 . These analytical results are in good agreement 
with the numerical results shown in Fig. |2]even for interme- 
diate A/kp lE9ll . 

Using the fermion Green function G(K), we can show that 
the fermion momentum distribution n(k) is isotropic and can 
be expressed as n(k) = (1/4) - i%/E%) & As shown 
in Fig. [3] with increased SOC, the distribution broadens, 
which indicates a BCS-BEC crossover. The new feature here 
is that the distribution generally displays nonmonotonic be- 
havior. The peak in the distribution is just located at k = A. 

The pair wave functions ^^(k) = (tAwA-iur') can ^ e 
evaluated as TT (k) = -(/A/4)e ie " a/E« and n (k) = 
-(A/4) 2« 1/^2 El, where e Wk = (k x + ik y )/\k\. There- 
fore, the superfluid state exhibits both singlet and triplet pair- 
ings for nonzero SOC. The numerical results for the ratio 
\<f>ll(k)\/\(f)ii(k)\ displayed in Figj3]show that the triplet pairing 
spreads to wider momentum regime with increased SOC. Ac- 
cording to the general formula for the condensate number of 
fermion pairs Q, No = 5 f J d 2 rd 2 r'\(tf/ CT (r)i// al (r')}\ 2 , 
the condensate density reads «o = Zktl^TlC 5 )! 2 + I0TtO^)I 2 ]- 



FIG. 3: (Color-online)(a),(b)&(c) The momentum distribution n(k) 
and the ratio R(k) = \(p^(k)\/\<p^(k)\ for various values of A/k F and 
fB/fF = 0.01. (d) The condensate fraction 2No/N as a function of 
A/kp for various values of 6b/6f- 



The triplet pairing amplitude contributes, in contrast to the 
fermionic superfiuids with only singlet pairing Dill . For large 
SOC, we find analytically that 2N /N = 1 - O(^) -> 1 (see 
also Fig. [3]), which indicates the Bose-Einstein condensation 
of weakly interacting rashbons. 

In the presence of a trap potential V(r) = jOJ^r 2 , the chemi- 
cal potential becomes p(r) = pc,-V(r) and the density distribu- 
tion n(r) can be solved from the constraint N = 2n J rdrn(r) in 
the local density approximation. As shown in Fig. |4j the atom 
cloud shrinks with increased SOC, which can be viewed as a 
preliminary experimental signal of the BCS-BEC crossover. 
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FIG. 4: (Color-online) The density profile n{r) (divided by nj = 
ep/7r) in presence of a trap potential for various values of A/k F . 
The Fermi energy = k\j2 in trapped system is defined as 6f = 
^[NTiojl (U, and the Thomas-Fermi radius reads Rj = -\/2ep/a>±. 



BKT transition temperature — At finite temperature in 2D 
we should rewrite the complex ordering field <t>(x) in terms of 
its modulus A(.y) and phase 0(x), i.e., O(x) = A(x) exp[i8(x)]. 
Since the random fluctuations of the phase 9(x) forbid long- 
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FIG. 5: (Color-online) (a) The superfluid density p s at T = (divided 
by ri) as a function of A/k 7 . The dashed lines represent the results of 
2m I 'iuq calculated from Eq. I0. (b) The BKT transition temperature 
as a function of A/k F . The dashed line represents the rashbon limit 
and the dash-dotted line is the mean-field result. 



range order in 2D, we have (<&(x)} = but (A(x)} + at 
T + 0. However, Berezinskii 13 33 an d Kosterlitz and Thou- 



less 13-411 showed that below a critical temperature Tbkt, there 
exist bound vortex-antivortex pairs and quasi-long-range or- 
der remains. 

To determine the BKT transition temperature, we de- 
rive an effective action for the U(l) phase field 0(x). 
To this end we make a gauge transformation tf/{x) = 
exp [i8(x)/2]x(x) |[l6i 171. Then we arrive at the ex- 



pression Z = / ADAD6 exp { - J?U e g [A(x), d0(x)}), where 
the effective action pU ef[ [A(x),d0(x)] = U' 1 J dxA 2 (x) - 
^TrlnS 1 [A(x), d6(x)] now depends on the modulus-phase 
variables. The Green function of the initial (charged) fermions 
takes a new form S" 1 [A(x), d0(x)] = @-\A(x)} - Z[d0(x)]. 
Here @~ l [A(x)] = G~ l [A(x), A(x)] is the green function of the 
neutral fermion, and E[<90] = t 3 [/3 t 6»/2 + (V#) 2 /8]-/[;V 2 6)/4+ 
NO ■ V/2] + (A/2)[T 3 cr x d y 6 - Icr y d x S\, where T t (i = 1 , 2, 3) are 
the Pauli matrices in the Nambu-Gor'kov space. 

Since the low-energy dynamics for A + is gov- 
erned by long-wavelength fluctuations of 6(x), we neglect 
the amplitude fluctuations and treat A as its saddle point 
value 1 3 171. Then the effective action can be decomposed 
as <Z/ eff [A(jc), d6(x)] 1/kiJA, d0(x)] + <Z/ pot (A). The potential 
part reads <Z/ pot /V = A 2 /U + £ k tf k -<W(E+)-'W(E h )] where 
= E/2 + Tln(l+e^ E ). The kinetic part can be obtained 
by the derivative expansion ptl^AA, d6{x)} = 27=1 ^Tr(^E)". 

Keeping only lowest-order derivatives of 9{x), we find 
that the kinetic term Hy m coincides with the classical spin 
XY-model, which has the continuum Hamiltonian Hxy - 
jST fd 2 r [V#(r)] 2 where the phase stiffness J = ^ and 
p s is the superfluid density 13511 . The superfluid density in 
our model can be evaluated as p s = n - p\ — p2, where 
Pi = (A/87T)Z a =±rdka({ a k + A 2 /&)[1 - 2f{E a k )]/E a k and 
p 2 = -(l/47r)2 ff=± 'kdk(k + aA) 2 f'(E a k ) The BKT 

transition temperature is determined by 7bkt = jSF 133-3^1. 



hand, for large ee and/or SOC, p s can be well approximated 
by its zero-temperature value for T ~ Tekt- We are interested 
in the case with small e B and large SOC. For large SOC, using 
the fact A <K \ju\, we find analytically that 112711 



Ps (T « T A ) - —n, J(T « T A ) - " 



— , (6) 



where n B = n/2 and m B is given by Eq. ©. Therefore, the 
phase stiffness J~ naturally recovers that for a Bose (rashbon) 
gas at large SOC. The BKT transition temperature and the 
phase stiffness jump AJ' reaches the rashbon limit Tbkt = 
mi B l(2m B ) = (2m/mB)ep/8 and AJ' = «b/»ib- To verify 
above analytical results, we show the numerical results for 
p s (T = 0) and Tbkt in Fig. [5] Even for weak attraction, a vis- 
ible pseudogap phase appears in the window r BKT < T < T A 
for A ~ k F . The SOC therefore provides a new mechanism for 
pseudogap formation in 2D fermionic systems. 

Finally, we point out a surprising result, p s < n at T — 0, 
which is in contrast to the result p. = n for fermionic superflu- 
ids in the absence of SOC HHQ. Actually, at T — 0, the su- 
perfluid density reads p s = n -px, where the /l-dependent term 
Pa — P\(T = 0) is always positive and is generally an increas- 
ing function of A, Therefore, the superfluid density shown in 
Fig. [3] has entirely different behavior in contrast to the con- 
densate density shown in Fig. [5] It is generally suppressed by 
the SOC effect. The exact two-body solution provides a very 
transparent explanation to this suppression. At large SOC, the 
effective mass m B > 2m is an increasing function of SOC and 
causes the suppression of the superfluid density by a factor 
2m/niB- Our argument also applies to the suppression of the 
radial (x — y plane) superfluid density pf for the 3D case 11211 . 
where the radial effective mass is larger than 2m ifioll . 
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Note Added — After finishing this Letter, we note that sim- 
ilar results of the condensate density Il2l 13811 and the super- 
fluid density Il2ll in spin-orbit coupled Fermi gases are also 
reported. 



For sufficiently small e B and SOC, A is correspondingly 
small and Tbkt recovers the mean-field result T A . On the other 
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Appendix: In this supplementary material, we present the derivation details of some results in the main text. 



(A) Two-Body Problem: Binding Energy and Effective Mass 

Using the free fermion propagators g±(K), T~ l (Q) can be expressed as 

1 1 



r-\Q) = - - - J] Tr [g + (K + Q)a y g-(K)cr y ] 

L K 



(7) 



Completing the Matsubara frequency sum, we obtain Eq. (2) of the text. For the two-body problem, we discard the Fermi-Dirac 
distribution function and define the solution for r _1 (o>, q) = as E q = -{to + 2p). The two-body equation becomes 



I 



2£kq 



k 2 + 



■ 4A 2 k 2 



4A*k 2 q 2 sin- tp 



= o. 



(8) 



Here <p is the angle between k and q, and £kq - E q + ek+q/2 + fk-q/2 = E q + k 2 + q 2 14. 

For zero center-of-mass momentum q, the above equation reduces to J Q kdk[2(k~ + 6b) - Za=±(* + 2a/lA: + Z?b) ] = 0. 
The integral can be carried out directly. The easiest way is to use the trick k 2 + 2Ak = (k ± A) 2 - A 2 . Since the integrals are 
logarithmically divergent, we can convert the integration variables to k + A. Finally we find that it becomes 



f 

Jo 



dz 



1 



■Eb) Jo k 2 + 



dk 



E B -A 2 



(9) 



\Z + €b Z + 

Using the condition E B > A 2 we then obtain Eq. (3) of the text. 

For nonzero center-of-mass momentum q, we write E q E B - q 2 /(2m B ) for small q 2 and expand Eq. ([8]) to the order 0(q 2 ), 
then we obtain 



2m\ r 
mb) Jo 



kdk- 



{k z + EbY + 4A z k 



2/,2 



[(k 2 + Eb) 2 - 4A 2 k 2 f 



f 

Jo 



kdk- 



8A 4 * 



4 1,2 



(k 2 + E B ) [(k 2 + E B ) 2 - 4A 2 k 2 f 



Defining k — Eb/A 2 , this equation becomes 



_ 2m _ r 03 8x r r 

m B Jo (x + k)[(x + k) 2 - 4x] 2 [Jo 



dx 



(x + k) 2 + 4x 



[(x + k) - 4xY 



Completing the integrals analytically, we obtain Eq. (4) of the text. 



(B) Derivation of the Ground-State Energy 



(10) 



(11) 



In the mean-field approximation, the ground-state energy can be expressed as 

\2 



U 2 P ^Y 

where the inverse fermion Green function reads 

I io)„ - ^k + hcr z - A(k v cr x - k x cr y ) 

y (ico„,k) - 



;cr v A 



-io-yA ia>„ + ^k - ho~~ - MJc y cr x + k x cr y ) 

Using the formula for block matrix, we first work out the determinant and obtain 

detfT^Wn, k) = [{icj n f + h 2 A 2 k 2 - A 2 f - 4h 2 (ito n ) 2 - 4A 2 k 2 - h 2 ) . 



(12) 



(13) 



(14) 



Then completing the Matsubara frequency sum and taking T — we obtain O = A 2 /U + (1/2) £ k (2& - E+ - E~) where the 
term 2k is added to recover the correct ground state energy for the normal state (A = 0). The quasiparticle dispersions are 
given by the positive roots of the equation det^ -1 = 0, i.e., 



^ + A 2 + A 2 k 2 + h 2 + 2 yjtl(A 2 k 2 + h 2 ) + h 2 A 2 



1/2 



(15) 
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For h — 0, they reduces to E^ — ± Ak) 2 + A 2 . At finite temperature, the thermodynamic potential reads Q. = 1/ po t/V = 
A 2 /t/ + - W(£+) - where *W(£) = E/2 + T ln(l + e^ £ ). 

For r = and h — 0, the ground-state energy can be expressed in terms of E B as Q = (A 2 /4n) Yia=± J kdk[(2ek + 2a Ak + 
£b ) 1 - (Ef + £?)]. Since the integrals are convergent, we can use the trick k 2 + 2Ak = (k± A) 2 - A 2 and convert the integration 
variables to k ± A. After a straightforward calculation, we obtain 

A 2 2A A 

Q. = Q 2 d(A, j u,£ , b)+ -: arctan - + Q^. (16) 

4 * Je^A 2 ^JEl^ r A 2 

Noticing the fact that Eg satisfies Eq. (3) of the text, we obtain Q = Q2d(A,//, e B ) + Qa- 



(C) Solution of the Gap and Number Equations at Large SOC 

The original forms of the gap and number equations at T — are 

4 + 



I = Iv _L + _L , n = y i-JL-JL 

U 2^\2E+ 2E-)' 2Et 2EZ 



2 K 2E k 



(17) 



For large SOC, we expect p < and A <K \fi\. Therefore, we can expand the equations in powers of A/|/i| and keep only the 
leading order terms. The gap equation becomes 



f 

Jo 



kdk 



1 



k 2 + e B ^ k 2 + 2a Ak - 2u 

D a=± ^ I 



0. 



We obtain // = -E B /2. Substituting this into the number equation, we obtain 



1 



" n 2ttZjJ (k 2 +2aAk + E B ) 2 



_ a2 r 

* Jo 



kdk 



(k 2 + E B ) 2 +4A 2 k 2 
[(k 2 + E B ) 2 - 4A 2 k 2 ] 2 



(18) 



(19) 



We notice that the integral also appears in Eq. ( TlOt . Completing the integral analytically, we obtain A = ~J2E B e F £(i<) where f(/c) 
is defined in the text. 



(D) The Fermion Green Function and Related Quantities 

The explicit form of the fermion Green function @(ia> n , k) can be evaluated using the formula for block matrix. For h = 0, we 
find that the matrix elements (in the Nambu-Gor'kov space) can be expressed as 

ky(T X - k X (Ty ky(T X + k X (Ty 

Q\\ = &\\ + ; Sli, @22 - &22 + ; '-^22, 



&i2 = -ia-y 



kycr x + k x cr x 

J{\2 + — ; -Si2 



@2l = i<T y 



k v cr x - k x cr v 



(20) 



Here Jlij and Sy take the forms 



a=± 



2^(icj„) 2 -(E?) 2 ' 



a\2 ' 



^22 



1 TTH 11 

= 2 hjM^n 



icj n - 



f-{Eir 



-^21 - -^12, 



(21) 



and 



#12 



~ 2^ ± a {ioj n ) 2 -{El) 2, 
1 y A 



22 2 ^ ± nu n ) 2 -{Eir 



21 - --D12- 



(22) 
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Using the matrix elements of the Green function, we can calculate various quantities. First, the momentum distribution can 
be evaluated as 



«(k) = <(A kT tA kT > = <<AkiiAm> = ^^n(/w„,ky w » 0+ 



Second, the singlet and triplet pairing amplitudes can be expressed as 



(23) 



<M k ) = <<AkT^-m> = ^ ^^2i(jw„,k), <* iT (k) = <(Aki*A-kT> = ^^2i('w„,k), 

£ v ~ 1 V~i + 1 V" 1 
n (k) = <lAkT<A-kT> = — — t g 2j S 2l0'w„, ^U( k ) = <^kl<A-kl) = — " £ 2j S 2l0'w„, k )- 



Therefore, we have the relations <pn(k) - — 0j/f(k) and ^ff(k) = — 0Ti(k). 

According to Leggett's definition J30ll . the condensate number of fermion pairs is given by 



(24) 



N 



\ Z f f ^VKMr)<Mr')>| 2 



(25) 



For systems with only singlet pairing, this recovers the usual result No - J J d 2 rd 2 r'\{tfff(r)i//i(r'))\ 2 . Converting this to the 
momentum space, we find that the condensate density no = No/V should be a sum of all absolute squares of the pairing 
amplitudes. The final result for T = is 



"0 = \ Z ku( k )| 2 + l^tWI 2 + I^TT( k )| 2 + l^li*)! 2 ] 
1 k 

= \lZ 



A 2 A 2 



(26) 



For large attraction and/or SOC, we expect A <K |//|. Using the number equation (JT7J and expanding all terms in powers of A/1/4 
we can show that 2Nq/N = 1 - <9(A 4 /M 4 )- Therefore, the condensate fraction approaches unity at large attraction and/or SOC. 



(E) Effective Action of the Phase Field 

To obtain the effective action for the phase field 6{x) to the order (V#) 2 , we notice that the available operators in E[cW] are 
Ei = T3(V0) 2 /8, T.2 = -TV8 ■ V/2 and Z3 = (A/2)[T3<T x d y 9 - lo- y d x 0\. According to the derivative expansion, we have carefully 
checked that there are four types of nonzero contributions: 



1/i ~ Tt(0Li), V 2 ~ Tr(^S 2 ^S 2 ), ^3 ~ Tr(£I 3 £I 3 ), V A ~ Tr(^Z 2 ^S 3 ). 



(27) 



Since the superfiuid state is isotropic, the phase stiffness should also be isotropic. We have carefully checked that all anisotropic 
terms vanish exactly. Completing the trace in the Nambu-Gor'kov and spin spaces, we finally obtain the following expressions 
for the four types of contributions: 



-z/ 4 



P n k J J 

J Z Z J K + S ' 1 + ^ 22 + ^ 22 + 2 ^ + 2S2 0] / d 2 r(W, 

-R Z Z T K + ^2 + ^l) f ^ r ( V0 ) 2 ' 
P n k W 

S Z Z T ^ii^ii - ^22^22 + 2J[ 2x S 2l ) f d 2 r(W) 2 . 

P n k 1 J 



(28) 
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Collecting all terms, the effective action is reduced to a spin XY-model Hamiltonian Hxy = \3 f d 2 r[V6(r)] 2 , where the phase 
stiffness J' is given by 

3 = \ Z Z [i (^ iie ' w0+ - ^ 2e - iw " 0+ ) + 1 K + s?, + ^ 2 + b| 2 + ue 2l + 2^) 



« k 



/8 

+ — (ytj, + + ) + "2" (-^11^11 - ^22^22 + 2^21^21) 

Completing the Matsubara frequency sum we then obtain the expression given in the text. 



(29) 



(F) Properties of the Superfluid Density 



First, setting A = 0, we find that p s = 0. Therefore p s vanishes exactly in the normal state, as expected. Second, for vanishing 
SOC, the expressions of p s and ST recover the well known form given in fl7ll . Here we will examine the behavior of p s for large 
SOC at T — 0. At zero temperature, the superfluid density reduces to 



p s =n-px, p,\ 



a r 

8^ Jo 



dk 



(30) 



Therefore, even at T — 0, the superfluid stiffness does not recover the result p s = n for ordinary fermionic superfluids. Let us 
show what happens at large A. In this case p -E%/2 and A <s \p\. Therefore, we can expand the expression in powers of A/\p\ 
and keep only the leading order terms. Doing so, we obtain (see Eq. (T% ) 



a 2 r 

n ^ —A 
8^r Jo 



kdk 



1 1 

+ 



(O 2 



a 2 r 

w Jo 



kdk- 



(k 2 + E B ) 2 + 4A 2 k 2 
[(k 2 + £ B ) 2 - 4A 2 k 2 ] 2 



(31) 



and 



Pa 



a 2 r 

* 8^ Jo 

a 2 r°° 

w Jo 



dk 



kdk- 



1 1 \ 1 


1 1 


£ rj 2 


k + ) 2 (O 2 


8^ 4 /t 2 




£ B ) [(^ 2 + E B )- 


--4^ 2 /t 2 ] 2 ' 



(32) 



Comparing the above results with Eq. ([Tol l, we find that p^/n = I - 2m /m B . Therefore, for large SOC, the superfluid density 
and the phase stiffness are reduced to 



2m 2m n n B 

Ps = — n, J= — — = — 



(33) 



where n% = n/2 is the density of rashbons. This means that, at large SOC, the phase stiffness self-consistently recovers that for 
a rashbon gas. 



